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We present a tetraquark interpretation of the charged bottomonium-like states Z^*^ (10610) and 
Z^{10650), observed by the Belle collaboration in the 7r^T(nS') (n = 1,2,3) and ■K'^hb{mP) 
(m =1,2) invariant mass spectra from the data taken near the peak of the T(5S'). In this framework, 
the underlying processes involve the production and decays of a vector tetraquark ¥(,(10890), e^e" ^• 
¥(,(10890) ^ [Zjf (10610)7r=F,Z±(10650)7r=F] followed by the decays [Z^* (10610), Zj±( 10650)] 
TT^T {nS),Tv^hi,{mP). Combining the contributions from the meson loops and an effective Hamilto- 
nian, we are able to reproduce the observed masses of the (10610) and Z,;^ (10650). The analysis 
presented here is in agreement with the Belle data and provides crucial tests of the tetraquark 
hypothesis. We also calculate the corresponding meson loop effects in the charm sector and find 
them dynamically suppressed. The charged charmonium-like states (3752) and Z^ (3882) can be 



searched for in the decays of the J 



PC 



tetraquark state F(4260) via ¥(4260) 



and ¥(4260) ^^ (3882)71^, with the subsequent decays (Z±(3752), Z±(3882)) 



Zt{3752)n'^ 
{J/^,h^)-K^. 



PACS numbers: 14.40Pq, 13.66Bc, 14.40.Rt 

Recently Belle [l| (updating a previous publication 
reported the measurement of the n^T{nS){n = 1,2,3) 
and Tr^ht{mP){m ~ 1,2) invariant mass spectra from 
the data taken near the peak of the T(55) resonance 
in the processes e+e" — >■ T(nS')7r+7r~ and e+e" — ?> 
hh{mP)TT~^ TT~ , in which two charged bottomonium-like 
states (10610) and (10650) are discovered. Here- 
after, these states will be abbreviated to Zi, and Z^, re- 
spectively. The masses and decay widths averaged over 
the five different final states are in MeV [l|: 



10607.2 ±2.0, 
18.4 ±2.4, 



m^,± = 10652.2 ± 1.5, 



11.5±2.2. 



The angular distribution analysis indicates that the 
quantum numbers of both Z^ and are I^{J^) = 
1+(1+). These states defy a standard bottomonium as- 
signment, as in the valence approximation they consist 
of four quarks bubd (and charge conjugates). 

Due to the proximity of the Zj, and masses with the 
BB* and B*B* thresholds it has been proposed that 
the former could be realized as S'-wave BB* and B*B* 
molecular states, respectively (3-[l0|. In this scenario, the 
heavy quark spin structure of the Z^, and is expected 
to mimic that of the corresponding meson pairs 
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where 0~ and 1~ denotes the para and ortho- states with 
negative parity, respectively. One anticipates the mass 
splitting to follow Am^t = "^z' ^"iZj, = ?ti_b» —?tib ~ 46 
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MeV, in neat agreement with the observed value Amz^ = 
(45±2.5) MeV Moreover, the structure in Eq. (P pre- 
dicts that Zt and Z^ should have the same decay width, 
which is approximately in agreement with the data. 

Despite these striking patterns, the fact that both Zi, 
and Zl lie nominally above their respective thresholds by 
about 2 MeV reveals a tension with the molecular inter- 
pretation. If consolidated by more precise experiments, 
this feature may become a serious problem in this ap- 
proach, as a one-pion exchange potential, which would 
produce such a bound state, does not support an S'-wave 
BB* resonance above threshold in an effective field the- 
ory [llj- Also, the measured total decay widths appear 
much too large compared to the naively expected ones 
for loosely bound states, and this suggests that both Zi, 
and Z'fj arc compact hadrons. In addition, the measured 
cross sections in question are too big to be interpreted in 
terms of the decays T{5S) — (T(nS'), /i;,(mP))7r+7r~. 

In this paper, we pursue a different ansatz in which the 
observed processes arise from the production and decays 
of a vector tetraquark Y;, (10890) |12Hl4j . having a (Breit- 
Wigner) resonant mass of [10888. 4t2;6(stat) ± 1.2(syst)] 
MeV and a width of [30.7t^;^(stat) ± 3.1(syst)] MeV 
[l^iIlBl- The mass and, in particular, the decay width 
of yf,(10890) differ from the Particle Data Group entries 
assigned to the T(55) 3. We propose that the states 
Zfc and ZIj seen in the decays of ¥{,(10890) are themselves 
charged tetraquark candidates having the flavor configu- 
ration [bu] \bd] (and charge conjugates) (see Refs. [13, [3l 
for earlier suggestions along these lines). Their neutral 
isospin counterparts with 1^ = have J^^ = 1"' and 
their masses were calculated in the effective Hamiltonian 
approach in [l^ . Ignoring the small isospin-breaking ef- 
fects [13, [l^ , ^b and Z'l^ have the same masses as those of 
their neutral counterparts. As shown below, these esti- 
mates yield a too large value for Amz^ compared to the 
Belle measurement. 

However, threshold effects and common decay chan- 



2 



ncls may play an important role beyond what can be 
described by the constituent quark model and its tran- 
scribed form adopted in [l^ to work out the spectro- 
scopic aspects for the tetraquark states. In particular, 
two hadronic states having the same quantum numbers 
may mix due to dynamical effects, leading to differences 
in their masses and decay widths. Typically, the result- 
ing mass shift is dominated by decays to the common 
states and reflects the partial widths to these states. A 
case in point here is the mass difference between the 
and D*^ , which is dominated by such common decay chan- 
nels, and is usually calculated by the meson loops, as dic- 
tated by the optical theorem [20] • Following essentially 
the same line of argument, we quantify this effect for 
the two charged-bottomonium-like states and Z'^^^. 
We recalculate the masses of Zi, and states by taking 
into account the meson loop contributions involving the 
Zweig-allowed two-body intermediate states BB*, B*B*, 
hh{mP)'K, T(nS)'K and rj^p. Theoretical estimates pre- 
sented here account for the observed masses; in particu- 
lar, the precisely measured mass difference Amz^ is re- 
produced in terms of the partial decay widths of Zh and 
Z'fj. This can be tested in future when the partial de- 
cay widths are measured precisely. In our approach, the 
mass eigenstates Zi, and are rotated with respect to 
the tetraquark spin states Zf, and Z^, and we determine 
this mixing angle. 

We also work out the corresponding meson loop effects 
for the charged charmonium-like states Z^ and Z'^^ , with 
each one of them belonging to an isotriplet. The masses 
of the electrically neutral states have been calculated 
and are predicted to have typical values m{Zc) = 3752 
MeV and m{Z'J = 3882 MeV Ignoring the smaU 

isospin-breaking effects, these estimates apply for the 
charged counterparts Z^ and Z'^ as well. We find 
that the meson-loop effects are in this case dynami- 
cally suppressed, as detailed below. However, we do 
expect that the production and decays of the Z^ and 
Z'^ will essentially mimic the patterns seen for their 
bottomonium counterparts Z^ and Z'^^ . In particular, 
Z^ and Z'^^, which are not measured so far, can be 
searched for in the decays of the neutral J^^ = 1^^ 
state y(4260), via the processes r(4260) Z±(3752)7rT 
and r(4260) Z±(3882)7rT, with the subsequent de- 
cays (Z±(3752),Z±(3882)) ^ (J/V-, /ic)7r=^- 

We start with the classification of the Zi, and 
tetraquark states in terms of the spin and orbital an- 
gular momentum of the constituent diquark and antidi- 
quark. The concept of diquark was introduced by Gell- 
mann in his epochal paper on quarks [2l| and since then 
has been widely discussed in the literature (for reviews 
on diquarks, see Refs. [l^, H^). A diquark has positive 
parity and may be a scalar (spin-0, or "good" diquark) or 
an axial- vector (spin-1, or "bad" diquark) |24| - [26| and is 
assumed to be a color antitriplet 3c. The states Zi, and 
arise from the production and decays of 1^6(10890), iden- 
tified with a linear combination of the two tetraquark 
states Y|;,g] = [fog] [65] {q = u,d) having the spin and or- 



bital momentum quantum numbers: S[i,q] = 0, S'jjg] = 0, 

^[ba\M\ = 0> L[bq][bq] = 1: ^'^^'^ t^e total Spin J[f,,][t^] = 

1 [12[. We shall iDe using a non-relativistic notation to 
characterize the tetraquark states •S'pqj; J[bq][bq]}-: 

which a matrix representation of the interpolating op- 
erators is used in terms of the 2x2 Pauli matrices ai 
{i = 1,2,3): 0[QQ] = Q^<72Q/V2, 1[qq] = Q'^(J2(j'Q/^ 
and Oqq = QQ/\/2, Iqq = Qa'-Q/y^, Q being any 
quark. The two tetraquark spin states Zt and Z^ are 
represented as 

\Zb) = {0[bq] ® - 1[M ® Olbq])/y^, 

\Z'b) =l[bq]®llbq]- (2) 

Performing a Fierz transformation, the fiavor and spin 
content in the bb ® qq and bq (8) qb product space can be 
made explicit: 

\Zb) = (-1,1 ® %q + ® l-,)/V2 = l-,^ ® l;,, 

\Zb) = (Ifel ® %q + ® lqq)/V2 

= i^bq®%l + %®\l)l^- (3) 

Eq. ([3]) shows that the and Z^ have similar coupling 
strengths with different final states. The labels Qbq and 
Ibq in Eq. ([3]) can be viewed as B and i?*, respectively. 
It follows that Zf, couples to B*B* state while Z^ couples 
to BB* . We stress that this identification is in contrast 
with the molecular interpretation, in which the transition 
Z^ — > BB* + h.c. is forbidden by the spin symmetry since 
Z^ is assumed to be essentially a B*B* molecule This 
difference can be tested in the future and is of great im- 
portance in order to distinguish between the tetraquark 
and the hadronic molecule interpretations. 

In the effective Hamiltonian approach, the 2x2 mass 
matrix for the 5'- wave l"*" tetraquarks M is given by [l^ 

where I is a 2 x 2 unit matrix, a = A/2 + (K^q)^ — K^q and 
6 = Kqq — K,,. lu the abovc (^65)3 accounts for the spin- 
spin interaction between the quarks inside the diquark 
and antidiquark, Kgg and Kjf, are the couplings accounting 
for the interaction between the the quarks in the diquark 
to the antiquarks in the antidiquark, and A is the mass 
difference between the spin-1 and spin-0 diquarks. Using 
the default values of the parameters [H, [iS] (in units of 
MeV) 

Kqq = 79.5, KijI = 9, Kbq = 5.75, (^65)3 = 6, (5) 

yields the diquark mass m^bq] — 5200 MeV (from the 
yf,(10890) mass). The value of A is uncertain, with A ~ 
200 MeV for the light quarks [1^. Reducing its value 
drastically for the c and b quarks will reduce the level 
spacing of the corresponding tetraquark states for which 
the experimental evidence is rather sparse. Due to the 
lack of data, we adopt an admittedly somewhat arbitrary 



range A = (120±30) McV for our numerical calculations. 
These parameters yield the following values for the two 
charged tetraquark masses and the mass difference 



mz, = (10443li^^)MeV 



(10628l^^)MeV, 

-.+21 



= (185+f^)MeV. (6) 



We note that the prediction for Amz^ given above is 
much larger than the experimental data, and there is no 
easy-fix for this mismatch at hand in terms of the pa- 
rameters in the effective Hamiltonian. Since this Hamil- 
tonian [l^ adequately describes the mass spectrum of 
the exotic states discovered in the charm sector, we con- 
tinue to use this as our starting point and argue that 
additional dynamical contributions to the mass matrix 
arise from the meson loops. 

With this premise, the renormalized masses can be ob- 
tained by computing the two-point functions. At the one- 
loop level, the self-energy corrections to the unperturbed 
propagator Yi{jP')gf^^, depicted in Fig. [U are written as 



Taking the /i;,7r state as an example, we find 



(47r)2 



dxsA 



(7) 



(8) 



where A = x'^s — xs + xrn^ + (1 — x)mf^^ — ie, and the 
coupling constants appearing above are defined through 
the hadronic interaction 

C = e^^^^g^^o^^^d'^zj^'^^d'^h^.TT + h.c. (9) 

In deriving S(s), the MS scheme in the unitarity gauge 
has been used to remove the UV divergence with the 
scale /i ^ m ^(,) . We recall that the real parts of 

E(s) contribute to the mass matrix, while the imag- 
inary parts of E(s) are related to the decay widths 
of Zb and Z^. In particular, the transitions Zb 
{T{nS)'K,hb{mP)7r,rib{nS)p) — > Z^ contribute to the off- 
diagonal terms in the 2x2 mass matrix and provide sig- 
nificant effects on the mixing of the two tctraquark-spin 
eigenstates. 

The meson loop corrections due to the different 
hadronic channels can be viewed as renormalizing the 
"bare" mass for the states predicted in the constituent 
tetraquark model. We are interested in the relative mass 
shifts which are reflected by the genuine part of the loop 



contributions ReEo 



These can be obtained b" 



1 DV 



subtraction procedure at a suitable mass scale sq |27| : 



Rel]ge„(s) = Rel](s) - ReI](so). 



(10) 



Setting the scale sq needs a prescription. It is reason- 
able to choose So as the mass squared of the lowest lying 
bound state for a given quark flavor, which we take as the 
J^'-^ = 0++ scalar tetraquark state. A different choice. 





T{nS), hij{mP),r]f,{nS) 

■K,p 



FIG. 1. Two-body meson-loop corrections to the function 
E(s) defined in Eq. ((8]). The intermediate states B^*^ B* con- 
tribute only to the diagonal terms in the mass matrix while 
T{nS)-K, hb{mP)-K and rji,{nS)p contribute to both the diag- 
onal and non-diagonal terms. 



namely sq = ^^fbqY ^'^^^ slightly modify our results and 
the effects caused by the ambiguity in sq will be incor- 
porated in estimating the systematic uncertainties in our 
approach. 

Including the loop corrections, we now have the follow- 
ing structure for the 2x2 mass matrix 



r. 



ZL 



(11) 



where i runs over the two-body channels shown in Fig. [T] 
the coefficients Ci(s, sq) are defined as 



R65]ge„(s) 

' 2 Iml](s) 



(12) 



in which s is taken as the physical mass squared from 



the data and Rel]gen(s) is determined as in Eq. (fTOt . 
The sign in the T{nS)TT contributions to the off-diagonal 
terms is reversed due to the spin symmetry as shown in 
Eq. In the case of open bottom mesons, the BB* loop 
impacts on M22 while B*B* modifies Mn. Note, that 
via the optical theorem the imaginary parts are directly 
related to the decay widths, and our parametrization in 
Eq. pT|) makes this manifest. 

Choosing the subtraction point as sq ~ [(10.385 ± 
0.05)GeV]2, which corresponds to the mass of the lowest 
(O"*""*") tetraquark state with a hidden bb quark content, 
we estimate the following values for the coefficients Ci 
(ignoring errors on the smaller qs): 



Chh(2P)7r 










-1.1 


3± 1 


-1.1 



For the analysis of T{nS)'!T contribution, the Lagrangian 
C = g^.^ioV^Z^li: with V = Yb, T{nS) gives 



ffz^'>T(«S)7r52(')T(„S)7r 



A 



(47r)2 



dx 



-log (4) 



A 



(13) 
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FIG. 2. Dependence of the coefficients a on tfie subtrac- 
tion scale So calculated with the Lagrangian in Eq. ((TSj. 
The frames (a) to (g) correspond to the intermediate states 
i = T(lS)7r,T(2S')7r,T(3S)7r,/ii,(lP)7r,/i6(2P)7r,r;6p,BB^, 
respectively. 



with A = x^s — xs + xm^ + (1 — a;)TOx(„5) — 
relevant coefficients are predicted as 

CT(lS)7r = -0.01, Cx(2S)7r = -0-1, 
CT(3S)7r = -1.3. 



.(') 



The expression for contributions from the ' — > i^bP 
channel is similarly obtained by replacing the vector 
TinS) by p and the pseudoscalar tt by yyf,. 

Instead of the Lagrangian specified above, using the 



Lagrangian with the derivative of the pion field which is 
inspired by the chiral symmetry 



(15) 



we have 



S(,s) = - 



Z^''T(nS)7r%<')T(«S)7r 



(47r)2 



dx 



A3 



4 log 



A 



- 5 



A2 



-8m 



T(nS) 

As 



3(4™^(„5)+5(2-3x2))log(4) 



- I2sx^ - 7s 

[{8x' - 4)m|(„5) - {x' ~ l)'s] 

A 



+ [(8 - 12.T2)m|(„s) + ,s(l - x'f] log (-) 

P 

-s^(.^-ifiog(A)|. 



(16) 



This expression yields larger values for the coefficients 

CT(nS)7r: 



CT(lS)7r 



-0.7, Cx(2S)7r = 
CT(3S)7r 



-2.1, 
-6.5. 



(17) 



It should be noted that these numbers are much larger 
than the ones in Eq. (|14|) . due to the fact that the pion 
momentum coming from the derivative in Eq. (|15p is 

small in the ^ rest frame and thus suppresses the par- 
tial decay width and hence the denominator in the defi- 
nition of Ci as in Eq. (|12p . 

The dependence of these coefficients on the subtraction 
scale is shown in Fig. [51 where the Lagrangian in Eq. ([15]) 
has been used. The striking result is that the coefficient 
c/i6(2P)7r dominates by far all the others. The main rea- 



son for the dominance of the coefficient Ch, 



is that 



the limited phase space and the p-wave decay character 
of Z^''' — )■ hb{2P)Tr result in a comparably small value for 
the i mag inary part of E(s) compared to its real part. In 
Ref. [25, Belle collaboration has reported the measure- 
ments of the cross sections for e~^e~ T{nS)TT^Tr^ and 
e"'"e~ — > hi,{mP)TT~^ 7T~ near the peak of the T{5S) res- 
ie and the onance. Using the final state T(2S')7r+7r~ as normaliza- 
tion, they found that the ratios of the various cross sec- 
tions are typically all of order 1. Thus, the partial widths 
for the different final states listed are comparable, which 
(14) suggests a value of 0(1) MeV for each of them [l|. Thus, 
the domination of the hi,(2P)T: channel in the meson-loop 
corrections to the 1^ mass matrix is a consequence of 
this channel having the largest coefficient and the antic- 



ipated sizable partial decay width of Zj, 
This is worked out quantitatively later. 



(') 



hb{2P)TT. 
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FIG. 3. Dependence of the coefficients Ci for hadronic decay 
channels of the charmonium-hke Zc states on the subtraction 
scale So calculated with the Lagrangian in Eq. (|15fl . The 
panels (a) to (e) correspond to the intermediate states i = 
J/i^TT, tp{2S)n, hcTT, rjcp, DD* , respectively. 



For comparison, we have performed the same calcula- 
tion for the hidden- charm tetraquark states whose masses 
are calculated in the constituent diquark model by Ma- 
iani et al. [l^, predicting the masses of the two 1"' cc 
hidden tetraquark states as 



3.752 GeV, mz' = 3.882 GcV. 



(18) 



This yields a mass difference lS.mz^ = 130 MeV. Ignor- 
ing the isospin symmetry breaking effects, typically a few 
MeV, the above estimates hold also for the charged coun- 
terparts. Since the above masses are very close to the 
estimate of the mass of the lightest scalar J^"^ = 0++ 
tetraquark state, ms^ = 3.723 GeV [l^, the genuine 
meson-loop contributions, after subtraction, are expected 
to be small. We show the corresponding coefficients for 
various hadronic channels in Fig.[3l where in order to de- 
termine the imaginary part the mass for the higher 1^ 
state has been used as the physical mass. A striking dif- 
ference between the coefficients shown in Figs. 2 and 3 




r::,2P).(MeV) (a) 



r(Zi,^/i6(2P);r)(MeV) (b) 



FIG. 4. Constrained partial decay widths from the Zb and 
Z'f, masses measured by Belle. The left-hand panel shows 
the constraint on the partial decay widths of the tetraquark 
eigenstates Zi, and Z',,. The circular (green) contour is ob- 
tained by the mass difference Amz^ = 45±2.5 MeV, while the 
slanted vertical (blue) band results from the averaged mass 
(mz^ + mz'^)/2 = 10629.7 ± 2.5 MeV for the defauh values 
A = 120 MeV and Ch^(2P)7r = 45. In the right-hand panel, the 
corresponding constraints on Zt and Z'^ partial decay widths 
are depicted. The solid (black) region results from default 
values, while the extended (red) region is obtained by varying 
these two parameters, as stated in the text. 



is the absence of the coefficient Ch^(2P)-n in Fig- 3, as the 

decays Zc'' — > /ic (2P)7r are not allowed kinematically. 
Assuming that the partial decay widths in the various 
channels J/tjjTr,^p{2S)Tr, hc{lP)Tr,r]cP etc. are of order 1 
MeV, as in the decays of the ^^'^ , we anticipate that the 
corrections due to the meson loops in the Zc - Z[. mix- 
ing are also typically of the same order, namely order 
1 MeV, hence not significant. Thus, unlike the masses 
of the Zfc -Z'^ complex, the masses for the hidden-charm 
tetraquarks calculated in Ref. [l^ are not expected to be 
significantly modified by meson- loop effects. 

As already remarked, the 1"' relatives of the Z^, and 
Z'y^ states in the charm sector, Zc and Z^, in our model are 
expected to be produced in the decays of a 1^^ hidden- 
charm tetraquark. The state V(4260) fits the biU. The 
enhancement of the cross sections for — )■ J /tpTT'^-K~ 
and e+e^ hcTT^n^ seen by the CLEO collaboration at 
the center-of-mass energy around 4.26 GeV [2^ is very 
likely a signature of their existence. In order to confirm 
or negate this scanrio, we suggest our experimental col- 
leagues to scan over this mass region more precisely. 

Returning to the discussion of the mass difference of 
the 1"*" tetraquarks in the hidden bb sector, we note that 
it is approximately given as Amz^ ~ 2\/ a'^ + b'^, where 

c,(r^-rf)/2, 6' 



a' = a - c,^i7" - r;'')/2, b' = b- cVrf'rf^ and i 
denotes hi,{2P)Tr, as we keep only the dominant contri- 
bution. The corresponding mass eigenstates are 



sin^zj^), 



\Zb} = cos 9 z,\Zb} 
\Zl) = sin 9 z,\Zk)+ cos 9 zjZ',), 

^-'[b'/{a' + Amzj2)]. 



(19) 



with 9zi, = tan 

In Fig.|31 we show the constrained partial decay widths 
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from the masses observed by Belle. The left panel shows 
the constraints on the widths of the tetraquark mass 

eigenstates zj^''^ for the default values of A and Ch^(2P)-K- 

z z' 

In the spin-symmetry limit, and are equal. As 
seen in this panel, the resulting contours intersect at two 
regions, the lower one of which implies large symmetry 
breaking effects and hence is not entertained any further. 
The upper region in which the two couplings differ by ap- 
proximately 40% is further analyzed. In the right-hand 
panel, the corresponding constraints on the Zi, and 
partial decay widths are depicted. The black region de- 
notes the default values of A and Ci given above, while 
the extended (red) region is obtained by the variations of 
these two parameters in the ranges A = (120 ± 30)MeV 
and Ch^{2P)TT = 45j^i[Q. Based on this, we estimate 

Oz, = (-19ti?)°, 

TiZb ^ h{2P)7r) = (1.7ilJ)MeV, 

T{Zl /ih(2P)7r) = (2.5l? j)MeV. (20) 

We note that the mixing angle is small, implying that the 
mass eigenstates are close to their respective tetraquark 
spin states. From the partial widths given above, we 
extract the relative strength of the coupling constants 

In the Belle data [l|, the ratio ?'hi,(2P)7r is not mea- 
sured directly; what is reported is the ratio Oie**^' = 
SYfcZ^TT X gz'^i/igy^ZtTT x gz^i), which are products of the 
production and the corresponding decay amplitudes of 
the Zh and in the given final states. The updated 
value in [l[ is a/i^(2P)7r = l-6lo;4to'.6- analysis to 
estimate the relative amplitudes in all five final states re- 
ported in Table I in [l| is in progress in the tetraquark 
context. We anticipate that the couplings in the pro- 
duction amplitudes involving Z^ and arc similar, i.e., 
ISnz^Trl - IgYbZiirl and hence r,,^(2P)^ = ah(,(2P)7r, in 
agreement with the Belle data. 

Using the Lagrangian given in Eqs. (jH]) and (ITSI) . 
we have the following amplitudes for the decays Yi, — > 
T(n5')7r+7r~ and Yf, — > h{,{mP)'K^ 

iM{Yb T(7i5)7r+7r") = A„Zi + ■igYiZt^gztT{nS)7r 

en ■ 4(„5){(BW^f"^)-" + ax(„5).e^^-'-'"BW^;""'-") 

(™n - PT{nS)n+)iPTinS)7T+ ' "4(nS)) + (^T^ ^ TT^)}, 

iM{Yb hb{mP)Ti+Ti^) = A'^^z^ - WYtZt,-KgZtK-,T(-^iuai3 

+a.,(,„P).e^^'-'"-"BW^f "^)^"] + {n+ ^ tt")}, (22) 

with BW* (,) = \pf — m'^,,, + im„(>)T ,^(,)]~^. Belle mea- 

surements show that the a^s are roughly 1 within errors, 
while the phases are close to either or 180°, though 
the errors are rather large. It is worth pointing out that 
the momentum dependence arising from the interaction 



Lagrangian given in Eq. (|15p are not taken into account 
in the parametrization adopted by Belle. Although the 
relative strength of the amplitudes, namely 0^6**^', is not 
affected, the 7r^T(nS') and Tr^hb{mP) spectrum distri- 
butions will be modified. Anz^ and A'^^z^ refer to the 
non-Z^''' amplitudes in the indicated final states. 

The structure of Anz^ was worked out in the 
tetraquark picture in great detail in Refs. [ij, [l^ . 
As opposed to the amplitudes involved in typical di- 
pionic heavy Quarkonia transitions, such as T(45) — > 
T(lS')7r+7r~, which are modeled after the Zwcig- 
suppressed QCD multipole expansion [30| . the ampli- 
tudes for the decays ¥(,(10890) — ?> T(nS')7r+7r~ are not 
Zweig-forbidden, and hence they are significant. In ad- 
dition, they lead to a resonant structure in the tttt in- 
variant mass spectrum. This is most marked in the 
T(lS')7r+7r~ mode in the form of the resonances /o(980) 
and /2(1270). The measured dipionic invariant mass 
spectra by Belle in these final states is in conformity with 
the predictions [Tsl . Il4^ . On the other hand, the ampli- 
tudes A'^z^ in the decays Fb(10890) — >■ hb{mP)Tr^n^ are 
expected to be neither resonant nor numerically signif- 
icant. Only the transition n(10890) ^ /ib(lP)/o(980) 
is marginally allowed, heavily suppressed by the phase 
space and the P-wave decay character. The state /o(600) 
(or cr(600)) contributes, in principle. However, as this is a 
very broad resonance, the higher mass part is again sup- 
pressed by the phase space and hence the contribution 
of the /o(600) in the decay ^(10890) /i6(lP)/o(600) 
is both small and difficult to discern. This feature is 
also in accord with the Belle data Q. Finally, we note 
that the absence of any anomalous production of the 
states (T(n5')7r+7r~, /ib(mP)7r+7r~) in the decays of the 
bottomonium state T (11020) [l^ is anticipated in the 
tetraquark picture, as opposed to the hadronic molecu- 
lar interpretation of the Z^ and for which the decays 

T(11020) ^ zl'^^TT^ T{nS)7r+n~,hbimP)7r+7r- are 
expected to be enhanced by the larger phase space com- 
pared to the corresponding decays from the T[bS). 

In summary, we have presented a tetraquark inter- 
pretation of the two observed states (10610) and 
Z^=^(10650). Comb ining the effective diquark-antidiquark 
Hamiltonian with the meson-loop induced effects, we 
are able to account for the observed masses in terms 
of the decay widths for the dominant channel Z^'''^ — 
hb{2P)TT^ , obtaining a ratio for the relative decay am- 
plitudes in the decays -Z'^'''^ hb{mP)'K^ which agrees 
with the Belle data. Together with the resonant tttt struc- 
ture in the decay modes ¥(,(10890) — >■ T(n5)7r+7r^, first 
worked out in [isl . [l3 |. this Letter provides additional 
support to the tetraquark hypothesis involving the states 
yf,(10890), Zj=^(10610) and 2^=^(10650). Precise spectro- 
scopic measurements foreseen at the Super-B factories 
and at the LHC will provide definitive answers to several 
issues raised here and will help resolve the current and 
long-standing puzzles in the exotic bottomonium sector. 
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